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^ ! 1 Introduction 

(N 



Most young stars are surrounded by accretion disks which are likely pro- 
genitors of planetary systems. The study of the physical structure, chemical 
i__ , composition, and evolution of these disks is therefore clearly a topic of great 

' interest. The thermal emission of dust grains in disks has provided a wealth of 

If) \ information on the disks' properties. At near-infrared wavelengths, the ther- 

mal emission is dominated by a small amount (~ 10~^ Mq) of warm material 
. , close to the star, and the emission is often optically thick. These wavelengths 

' serve as excellent probes of the presence or absence of a disk. At the longer 

I (sub) millimeter wavelengths, the emission probes the cold material at larger 

radii (> 1 AU), where the emission is usually mostly optically thin. These 
wavelengths are good probes of the total mass of a disk. 

Dust, however, only makes up 1% of the total mass of the disk, if we 
adopt the standard interstellar gas-to-dust mass ratio. Ignoring the difficult 
to observe H2 molecule (but see, e.g., [2 El)) circumstellar disks are now 
commonly observed through lines of molecules such as CO, HCO+, HCN, 
CN, etcetera. Although it is found that the molecular abundances are often 
O I strongly depleted through freezing out onto dust grains and the lines 

^ ■ correspondingly weaker, molecular-line observations offer a powerful array 

^ I of diagnostics to study disks: from the relative intensities of transitions, the 

molecular excitation can be derived and the density and temperature of the 
gas constrained; from the line strengths, the molecular abundances are found 
I and the chemical composition inferred; and from the resolved line profiles, 

I the velocity field in the disk is directly probed. This contribution discusses 

the methods with which this wealth of information can be extracted from 
molecular lines, and presents examples of recent research that illustrates this. 



2 The Physical and Chemical Structure of Disks 

Within the framework of standard accretion disks 01 jSl , recent theoretical 
models specifically address their vertical structure [Hl[71- In a system where 
the central object dominates the mass and the luminosity, the disk is expected 
to assume a flared configuration, where the hydrostatic scale height increases 
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with radius [SHU]. This flared surface intercepts more steUar light, which raises 
the temperature of the surface and, through reprocessing, of the disk interior 
(e.g., [7]). Refinements to these models include the effects of dust sublimation 
and self-shadowing of the outer disk by inner, flared regions |1(J[I11[[T^I13| . 

The resulting descriptions of the density and temperature as function of 
radius and height then serve as input to calculate the abundances of var- 
ious molecules (e.g., [H HSl HEl El EB. 19, IQ, IL HJi; and Henning, 
this volume). Three processes affect the molecular gas- phase abundances (re- 
viewed in ^1): (1) freeze-out of molecules onto dust grains in the dense 
and cold midplane, reducing the gas-phase abundances; (2) evaporation of 
these ice mantles at intermediate vertical heights and near the star, where 
the temperatures exceed the evaporation values (~ 22 K for CO, 90 K for 
H2O) and where species return to the gas; and (3) dissociation and ionisa- 
tion by (inter-) stellar ultraviolet photons in the disk surface where there is 
insufficient shielding by dust, H2, and CO, and where a chemistry akin to a 
photon-dominated region is set up (see Sternberg, this volume). Radial inflow 
and vertical convection further influence the disk's chemical structure. 

3 Radiative Transfer and Molecular Excitation 

With a model in hand for the physical and chemical structure of the disk 
(density, temperature, systematic and turbulent velocity fields, and gas-phase 
abundances - all as functions of radius R and height z), the emergent emis- 
sion (or absorption) lines can be calculated. An iterative approach is usually 
chosen: Starting with an educated guess for the level populations throughout 
the disk (or at a sufficient number of grid points), the absorption and emis- 
sion coefficients are evaluated and the radiative transfer through the disk is 
solved (along a sufficient number of photon paths). This yields the intensity 
of the radiation field throughout the disk, with which the equations of statis- 
tical equilibrium can be solved and an update to the level populations found. 
This procedure is repeated until the solution converges. 

Several factors complicate this prescription in circumstellar disks. Disks 
have a large range in density and temperature, with molecular excitation 
ranging from thermodynamic equilibrium to sub-thermal. Molecular abun- 
dances can vary by several orders of magnitude within the disk. And opaci- 
ties, both in continuum and in lines, can be large. These factors require that 
the input model contains sufficient detail to cover the variations in excita- 
tion and abundance, and necessitate careful assessment of the convergence of 
the solution: large opacities can strongly slow down the convergence, because 
at each iteration step, changes in the population only propagate over short 
distances (corresponding to r ~ 1). 

Increased computer speed and memory has made realistic model calcu- 
lations tractable. In 1999 in a workshop at the Lorentz Centre in Leiden, 
several of the leading codes for molecular excitation and radiative transfer 
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were (successfully!) tested against each other^. The reliability of the obtained 
solutions depends on the quality of the molecular collision rates that are used; 
several projects are under way to provide accurate rates for astrophysically 
relevant species (see Roueff, this volume; and Schoier et al. in prep.^). 

Slow convergence of iterative methods in the presence of large opacities is 
not only an inconvenience for the impatient, but may lead to incorrect results. 
Convergence is often judged on the difference between subsequent iterations, 
with the implicit assumption that this will be small only when the solution 
is close to the correct value. The small differences between subsequent solu- 
tions due to opacity is easily misinterpreted as convergence. In the modeling 
of stellar atmospheres these problems are well understood, and have lead 
to the development of acceleration techniques such as Accelerated Lambda 
Iteration ,2^ (ALI; Lambda Iteration refers to iterative solution methods). 
Acceleration techniques have been included in Monte Carlo approaches to 
solve the radiative transfer (2^1 ^- This combines the reliability of accelerated 
convergence with flexibility in addressing arbitrary geometries. 



4 Examples 

Resolved Disk Emission of LkCa 15. Using the millimeter array of the Owens 
Valley Radio Observatory (OVRO), the emission of the 450 AU radius, 
0.18 M© disk around the young star LkCa 15 has been resolved [23 1^ . 
""^^CO J=2-l emission peaks at the source, with velocities consistent with 
Keplerian rotation of an inclined {i « 60°) disk around a 0.8 M© star. In 
stark contrast, the emission of CN 1-0 and HCN 1-0 is broken up in two 
peaks, symmetrically displaced from the star by ^ 2" w 300 AU. 

The appearance of HCN and CN can be reproduced by a simple radial 
'step function' for the abundance, with CN and HCN only present outside 
300 AU j29j . This can be explained by chemistry in a flared disk, where only 
outside ~ 300 AU densities are low enough to allow ultraviolet radiation to 
penetrate and set up a PDR with increased abundances of HCN and CN. 
Self-consistent calculations are currently under way (Kessler et al., in prep.). 

Two-dimensional Transport of Ultraviolet Radiation. To accurately account 
for the penetration of ultraviolet radiation from the central star into the 
flared disk, the full details of the two-dimensional transfer must be included 
|30|. Compared to one-dimcnsional calculations that only include vertical 
transport, the predominantly forward-scattering nature of the grains causes 
the ultraviolet field to penetrate deeper. A larger fraction of the disk is subject 
to photo-chemistry, with increased abundances of radicals like CN and C2H. 

^ For more information, see http://www.strw.leidenuniv.nl/~radtrans 
^Available through http://www.strw.leidenuniv.nl/~moldata and 

|http : / /basecol . obs-besancon. f r 1 

^ Available at http: //www. strw. leidenuniv.nl/~michiel/ratran.html] 
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An interesting application is the effect of the accretion-generated ultra- 
violet radiation field on the HCN and CN chemistry. Two-dimensional cal- 
culations show that the CN/HCN abundance ratio increases by a factor of 
10 if an ultraviolet excess similar to that of TW Hya is included on top of 
a 4000 K stellar blackbody [301 ■ The ratio of the CN 3-2 line emission over 
HCN 4-3 increases by a similar amount, albeit confined to ~ 1" region from 
the star at a typical distance of 140 pc for the object. Instruments like the 
Smithsonian Submillimeter Array (SMA) and the Atacama Large Millimeter 
Array (ALMA) will be able to trace such regions. 

Far-infrared CO Line Emission from Superheated Disk Surfaces. The object 
Elias 29 is a T Tauri star with a face-on disk of 0.012 Mq that is obscured 
behind several layers of interstellar clouds [HJj- It shows ^^CO 6-5 and ^^CO 
lines up to at least 20-19, indicating the presence of ^ 0.001 M0 of warm 
(> 200 K) gas [32] ■ This reservoir of warm gas may be explained by material 
heated in outflow shocks, but is also comparable to the amount of material 
expected in the superheated surface layer of the flared disk around Elias 29. 

Standard flared disk models [J] cannot produce sufficient far-infrared 
^^CO line emission: although the amount of warm gas is sufhcient, its density 
is too low to excite these lines [321 ■ standard flared-disk model, perfect 

mixing between gas and dust is assumed. If, however, the grains have settled 
toward the midplane, the superheated layer extends further into the denser 
disk interior. Calculations indicate that ^^CO far-infrared line strengths in 
this case can even exceed the observed values. Although more detailed model- 
ing is required, including the fate of small dust grains that will not decouple 
from the gas, this scenario suggests that the superheated layer is a viable 
source of far-infrared CO emission. 

Rotation and Infall in the Disk Around L1489 IRS. Submillimeter-continuum 
and millimeter aperture-synthesis observations show that the object L1489 IRS 
is surrounded by a flattened structure of ~ 2000 AU radius [331 ^.nd 
^ 0.02 M0 mass [331. While this mass is entirely typical for disks around 
T Tauri stars, its size exceeds by a factor of 2-3 the largest known gas disks. 
Through modeling of the HCO"'" 1-0 and 3-2 interferometer images, the ve- 
locity field in the disk can be obtained S3;- In addition to a Keplerian com- 
ponent, suggesting a central mass of 0.65 Mq, inward motions amounting to 
~ 10% of the total velocity vector are found. If such motions continue inward 
to the star (the interferometer data only sample scales > 700 AU), the disk 
lifetime is ~ 2 x 10'' yr. This is much shorter than the estimated duration 
of the embedded (few times 10^ yr) and T Tauri (few times 10^ yr) phases. 
It suggests that L1489 IRS represents a short-lived transitional stage, where 
the collapsing cloud core is settling onto a rotationally supported disk. 

The inward motions can be followed to within 0.1 AU from the star 
through ^^CO fundamental ro- vibrational absorption lines against the stellar 
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continuum at 4.7 /zm ,35^. The line profiles of the individual components of 
the P- and R- branches show prominent wings extending to +100 km s^^. 
Such velocities are expected to occur within 0.1 AU from the star, based on 
the velocity model derived from the HCO"'' measurements above. The pres- 
ence of absorption line originating from rotational levels of J— 6 and above 
indicate temperature exceeding 100-200 K, which also are only expected to 
be populated close to the star. Taken together, this shows that the inward 
motions are present from 2000 AU as traced by the HCO+ data to within 
0.1 AU probed by the CO absorption. 

There are several objections against the interpretation that the entire disk 
around L1489 IRS is contracting. The implied mass accretion rate exceeds 
observational limits by more than an order of magnitude. And the theo- 
retical models for accretion flow predict subsonic motions, while the observed 
speeds are supersonic. A possible solution is that only a thin, unstable sur- 
face layer participates in the inflow, while the bulk of the disk is Keplerian. 
Another solution posits the existence of two regions of inflow: one at large 
radii (probed by the HCO+ data) where the cloud's material spirals onto the 
outer disk, and one at small radii (probed by the CO line wings) where ma- 
terial streams from the disk onto the star. The fact that both regions can be 
fit with the same velocity model does not require them to be physically con- 
nected, because the motions occur in the gravitational potential of the same 
star. Further study of the L1489 IRS aims at distinguishing these scenarios 
(Hogerheijde et al., in preparation). 

5 Conclusion 

In this contribution I have shown that there now exist increasingly realistic 
two-dimensional and three-dimensional models of the physical and chemical 
structure of disks. Radiative transfer and molecular excitation tools have been 
developed that can calculate the emergent line spectrum on the basis of these 
disks models. On the observational side, current data are beginning to probe 
the structure of disks, although with no more than a few resolution elements 
across the objects. The first results show interesting structural differences 
between various probe molecules, which current chemical models can explain. 

In the next few years, increasing detail is to be expected in disk mod- 
els, as the increased speed of computers will bring more extensive modeling 
within reach. With the arrival of new observational facilities in the submil- 
hmeter (SMA, CARMA, ALMA) and the infrared (SIRTF, SOFIA, Herschel), 
molecular-line modeling will form an essential link in comparing models with 
data. By the end of this decade, our understanding of protoplanetary disks 
will likely have progressed significantly. 

MRH wishes to thank the organisers of the conference and the Leids 
Kerkhoven-Bosscha Fonds for their generous financial support. 
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